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Abstract 

Results on diffractive photoproduction of 4'{2S) mesons are presented using data col- 
lected between 1996 and 2000 with the HI detector at the HERA ep collider. The data 
correspond to an integrated luminosity of 77 pb^^. The energy dependence of the 
diffractive tp{2S) cross section is found to be similar to or possibly somewhat steeper than 
that for J/^ mesons. The dependences of the elastic and proton dissociative 
tp{2S) photoproduction cross sections on the squared momentum transfer t at the proton 
vertex are measured. The t-dependence of the elastic channel, parametrised as e**, yields 
^J2S) ^ ^ ^ Q^g^j compatible with that of the J/ip. For the proton 

dissociative channel the resuh b'^jf^^ = (0.59 ± 0.13 ± 0.12) GeV-^ is 2.3 standard de- 
viations smaller than that measured for the J/ip. With proper account of the individual 
wavefunctions theoretical predictions based on perturbative QCD are found to describe the 
measurements well. 
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1 Introduction 



We report on a measurement of diffractive photoproduction of ipi^S) mesons, 7p — > ^(25) +Y , 
where Y denotes either a proton or a proton dissociation system of mass My > nip. The data 
were taken with the HI Detector at the ep collider HERA. The dependences of the cross section 
on W^p, the photon-proton center of mass energy, and on t, the square of the four-momentum 
transfer at the proton vertex, are measured for the first time and are compared to those of J / 
production. The data cover the region 40 < W^p < 150 GeV. 

In recent years diffractive production of J/'ip mesons has been measured at HERA with increas- 
ing precision [|l]-§p. The cross section for the elastic process 'jp J/tjjp was found to rise 
steeply with W.yp. This was interpreted as a signature for a "hard" process and calculations in 
perturbative Quantum Chromodynamics (pQCD) were able to give a good description of the 
data In perturbative QCD, diffractive charmonium production in 7p scattering can be 

viewed in the proton rest frame as a sequence of several steps separated in time. An almost 
real photon is emitted from the incoming lepton and fluctuates into a cc pair. The cc pair subse- 
quently interacts with the proton via the exchange of two gluons (or a gluon ladder) in a colour- 
singlet state and then evolves into a real vector meson. In such models, with the cc fluctuation 
of the photon treated as a colour dipole, a number of distinct predictions are made for ^{2S) 
photoproduction [[7|-|TT|]. The cross section for ^{2S) should be suppressed with respect to that 
for J/ if), the dependence on W^p should be somewhat steeper than that of the J / if) photopro- 
duction cross section and the t-dependence of ^^(25) production should be similar or somewhat 
shallower than that of J /if) production. These predictions take into account the ^^(25*) meson 
wavefunction which is different from that of the J/^ meson in two respects. It has a larger 
expectation value for the radius than the ground state and it has a node (see e.g. [[T0|]). The 
cross section for elastic ^(25") production has been measured previously at HERA energies and 
the resulting ratio (T^(25)/aj/^ = 0.150 ± 0.027(stat) ± 0.018(syst) ± 0.011(fii?) [|l|] verifies 
the prediction of a suppression with respect to J /if) production. This and the other predictions 
mentioned above are addressed in this paper. 

In a geometric interpretation one expects the t-dependences, parametrised as oc e^*, to reflect 
the sizes of the interacting objects, similarly to the hadroproduction case. Vector mesons with 
large radii should have larger values of h than vector mesons with small radii. This has indeed 
been confirmed in photoproduction in the HERA energy range, where e.g. 6pO > was 
measured [|0],[T^]. Since the radius of the ^^(25*) meson is approximately a factor of two larger 
than that of the J meson one might naively expect a steeper t dependence. In QCD a different 
result is obtained. Due to cancellations in the contributions to the production amplitude from cc- 
quark dipoles with sizes above and below the node of the ^/^ (25) wavefunction, the t-dependence 
of elastic il){2S) production has been predicted to be slightly shallower than that of the J 
meson 

Calculations of the t-dependences for light and heavy vector meson production also exist [ [T3| ] 
using an additive quark model ansatz. These calculations are able to reproduce the measure- 
ments of the t-dependences at small \t\ for many elastic and proton dissociative processes. 

2 Data Analysis 
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2.1 Kinematics 



The kinematics of the process ep — * eipY are described by the following variables: the square 
of the ep center-of-mass energy s = {p+ kY; the negative four-momentum transfer squared at 
the lepton vertex Q"^ = —q^ = —{k — k'Y; the four-momentum transfer squared at the proton 
vertex t = (p — p'Y and the scaled energy transfer y = {p ■ q)/{p ■ k). The four- momenta k, 
k', p, p' and q refer to the incident and scattered lepton, the incoming and outgoing proton or 
excited state Y and the exchanged photon, respectively. The elasticity z of the meson production 
process is defined as z = (p ■ p^) /{p ■ q) where p^ denotes the four-momentum of the produced 
vector meson. In the proton rest frame, z describes the fractional photon energy transferred to 
the vector meson. It is related to My by ^ ~ 1 — (My —m?p — t) /W^p. For elastic events at low 
in which the proton stays intact, 2; ~ 1. 

In the photoproduction domain, i.e. at — > 0, the incoming lepton is scattered at small angles 
and below ~ 1 GeV^ it is not observed in the central detector. In the limit of photoproduction 
W^p = {p + qY is given by ys where y = {E — pz)^/{2Ee). Here E and pz denote the total 
energy and momentum component of the vector meson parallel to the proton beam direction|] 
of the vector meson and E^. is the energy of the incident beam lepton. For the measurement 
of the elasticity the relation 2; ~ (£' — pz)tp/'E{E — p^) is used where T.{E — p^) includes 
all measured final state particles. The variable t is approximated by the negative transverse 
momentum squared of the vector meson, i.e. t ~ —p't-s- 



2.2 Detector and Data Selection 



The HI Detector is described in detail elsewhere [|T^]. For this analysis the central and for- 
ward tracking detectors, consisting of a system of drift and proportional chambers with a polar 
angular coverage between 7° and 165°, are used for the detection of charged decay particles. 
In addition the main and backward calorimeters, covering the polar angular regions 4° — 153° 
and 153° — 177.5°, respectively, are used for lepton identification and for the determination 
of the event kinematics [|16|,|17]]. The instrumented iron return yoke of the solenoidal magnet 
(4° < 6* < 171°) which surrounds the central HI detector supplements the calorimeters in the 
identification of muons. Triggers based on lepton and track signatures are used to collect the 
events. 

The analysis presented here is based on data corresponding to an integrated luminosity of 
77 pb^^. The data were taken in the years 1996 to 2000. Until 1997 HERA was operated 
with positrons of energy 27.5 GeV and protons of 820 GeV, while afterwards the proton en- 
ergy was 920 GeV and both electrons and positrons were used. Events from il){2S) and J /if) 
production are selected using the direct decays into two leptons or e+e^ and the ^^(25*) 

is also reconstructed via the cascade into a J /if) and two charged pions with subsequent decay 
of the J/i}) into two leptons. The data selection resembles closely the procedure used in p,[T^]. 
Events with a scattered lepton candidate detected in the calorimeters with an energy deposit of 
more than 8 GeV are rejected. The accepted photoproduction event sample covers the range 
< 1 GeV^ with an average (Q^) ~ 0.055 GeV^ as determined from the Monte Carlo simu- 
lation. 



'The coordinate system of HI defines the positive z-axis to be in the direction of the proton beam. The polar 
angle 9 is then defined with respect to this axis. 
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Figure 1: Mass spectrum {\t\ > 0.07 GeV^) for the direct decay channel into leptons (/i^yU" and 
e+e") after the final selection. The insert shows the mass distribution restricted to the tpC^S) 
region. The curve shows the result of a fit combining Gaussian distributions for the J/ and 
ip{2S) signals, an exponential parametrisation of the radiative tail in the electron decay channel 
for the J/iJj and a linear background. 

Vector meson decays into lepton pairs are selected by requiring exactly two tracks in the central 
tracking chambers, each with a transverse momentum greater than 0.8 GeV in a polar angular 
rangeQof 20° < 9 < 160°. The tracks must define an event vertex in the ep interaction region. 
Both tracks are required to satisfy lepton identification requirements. The decay electrons are 
identified using the electromagnetic section of the calorimeters and energy loss dE/dx in the 
tracking chambers. Muons are identified in the instrumented iron return yoke or as minimum 
ionizing particles in the main calorimeter. Cosmic ray muon events are rejected using a track 
acoUinearity requirement. For the selection of ^'(25') events via the cascade decay, ip{2S) 
{J/ip — > £"'"£^)7r+7r^, events with exactly four tracks are selected. In addition to two lepton 
candidate tracks, which are identified as described above, exactly two further central tracks 
with opposite charge are required, each with a transverse momentum of at least 0.12 GeV. The 
invariant mass of the two leptons is restricted to 2.4 (2.9) < M« < 3.3 GeV for cascade decays 
with J/^/' — > e^e~ (J/ip yU+/i~), respectively. 

For the selection of the final diffractive (i.e. the sum of elastic and proton dissociative) event 
samples, further requirements are applied which ensure that the detector is essentially "empty" 
apart from the J/-ip and "ipi^S) decay products. These requirements and the method used to 
separate the elastic and proton dissociative contributions are discussed in detail in section 
Figure [1| shows the di-lepton mass distribution for the diffractive sample in the mass region of 

^For the data sample recorded in the year 2000, the polar angle of the decay electrons was required to be in the 
range 30° < 6* < 150°. 
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Figure 2: Mass spectra for the selected i£nn sample; a) and b) show the ii and iimr invari- 
ant mass distributions; c) shows the mass difference Mu-^t, — Mu when selecting the ii pairs 
reconstructed in the J/ if) mass window. 



the J and ^^(25*) mesons. The insert shows the same mass distribution restricted to the ^^(25*) 
region revealing a clear signal of 307 ±27 events above the non-resonant background. In decays 
with electrons the signals have a tail to lower masses due to the presence of radiative decays 
and radiation in the material of the detector. The non-resonant background is dominated by the 
process 77 in which the two photons couple to the beam lepton and to the proton, 

respectively. For the direct decays of tpC^S) mesons the non-resonant background fraction is 
large. It is smaller for the J/'ip^ events and for the cascade decays of ^^(25) mesons. 
Figure ^ and b show the £i and iinn mass distributions for the 4-prong sample. In Fig. |^ the 
mass difference Mut^t, — M« is shown. A signal of 278 events in the range iM^^r^ — Mu — 
0.59 GeV| < 0.06 GeV is seen with negligible background. 



2.3 Monte Carlo Models and Acceptances 

The acceptances and efficiencies for triggering, track reconstruction, event selection and lepton 
identification are calculated using Monte Carlo simulations. The J/ip and ^^(25*) samples are 



generated using the program DIFFVM [ [T8| ] and are passed through a detailed simulation of the 
detector response based on the GEANT program [ [TP| ] and the same reconstruction software as 
was used for the data. DIFFVM generates events according to the cross section dependences 
da/dt oc W^pC^^ for elastic photoproduction of charmonium. For production with proton disso- 
ciation d'^a/dtdMy oc W^pe^'^My is used. The parameters were chosen such that the main fea- 
tures of J/iIj and \Ij{2S) production are described: h = 4.8 GeV^^; b' = 1.6 GeV^^; /3 = _2.16 
and 4e = 0.96. Possible deviations from these parametrisations are taken into account in the 
systematic error analysis. The decay angular distributions of charmonium decaying directly into 
two leptons are simulated assuming s-channel helicity conservation. For background estimation 
the generators LPAIR [ PD| ] and GRAPE are used, which simulate the process 77 



The contributions of radiative decays to J/ijj or iIj{2S)^ e+e are estimated with the generator 
PHOTOS [O]. 



Detailed comparisons between the simulation and the data are facilitated by the large J/ip data 
sample in which the background fraction is small. Figure |3| shows the transverse momentum and 
polar angular distributions of the tracks in the direct decays of J/tp into muons and electrons 
both for data and simulation. There is good agreement. Figure ^ shows that the simulation 
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Figure 3: Distributions for the decays of J /if) mesons into muons (a,b) and electrons (c,d)- The 
polar angles (a,c) and the transverse momenta (b,d) of the decay leptons are shown for the data 
(points) and the DIFFVM Monte Carlo simulation (histograms). The Monte Carlo simulation 
is normalised to the number of events in the data. 



for the cascade decays of V'(2S') adequately describes the data. The simulation generates a 
flat TTTT-mass distribution in the available phase space, ignoring the dependence arising from 
the matrix element. To correct for this the simulation was reweighted using an event weight 
oc (M^^ — AmIY as determined in p3|]. This results in a good description of the Af^Tr spectrum 
(Fig. W)- 

The detector simulation has been checked extensively and separately for each data taking pe- 
riod and decay lepton type using independent data samples and was adjusted where necessary. 
Typical trigger efficiencies are 45% for the muon channel and 55% for the electron channel. 
The efficiency for identifying a muon (electron) is typically 75% (85%). Together with the 
geometric detector acceptance the overall efficiency varies between 5 and 10% with W^p. In 
the range studied here the variation of the overall efficiency with \t\ is small. The resolu- 
tion in t for the muon decay channel is typically 0.035 GeV^ and increases to 0.060 GeV^ at 
\t\ = 1.2 GeV^. For the electron decay channel the resolution is roughly 15% worse than this 
due to bremsstrahlung. 



2.4 Separation of Elastic and Proton Dissociative Processes 

For the analysis of the t-dependence it is necessary to separate elastic events from those with 
dissociation of the proton into a small mass system. To tag events with proton dissociation the 
forward section of the calorimeter (9 < 10°), the Proton Remnant Tagger (0.06° < 9 < 0.26°) 
and the Forward Muon Detector (3° < < 17°) are used. With these forward detectors, disso- 
ciated proton states with masses My > 1.6 GeV can be tagged [^. In the following, events 
with (without) signals in the forward detectors are called 'tagged' ('untagged'), respectively. 
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Figure 4: Distributions for the cascade decays of ^/'(2S') mesons: a) the transverse momentum 
of the decay leptons and b) the decay pions; c) the polar angle of the decay leptons and d) the 
decay pions; e) the two-pion invariant mass and f) the photon proton center-of-mass energy 
Wjp. The points show the data. The histograms show the DIFFVM Monte Carlo simulation 
normalised to the number of events in the data after event reweighting (see text). 



In the 'untagged' sample no additional tracks other than those from the J/ip or ip{2S) decay 
products are allowed. In the 'tagged' sample at most one additional track measured in the 
forward tracking chambers in the polar angular region below 10° is allowed and an elasticity 
z > 0.95 is required in order to ensure that the event sample is diffractive. 

The elastic and proton dissociative contributions to the diffractive sample are extracted from 
the untagged and tagged event samples respectively, taking the background admixtures from 
proton dissociative events in the untagged sample (~ 15%) and elastic events in the tagged 
sample (~ 10%) into account. The admixtures are due to the limited angular coverage of the 
forward detectors, their inefficiencies and noise fluctuations. 

In Table [T] the estimated relative signal and background contributions are listed for the tagged 
and untagged samples of the different decay channels. The non-resonant background in the 
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i,{2S) 




^(2^) - 










untagged 


tagged 


untagged 


tagged 


untagged 


tagged 


Elastic 


39% 


6% 


85% 


10% 


80% 


9% 


Proton Dissociative 


6% 


54% 


15% 


90% 


14% 


85% 


Non-resonant Bg. 


55% 


40% 






6% 


6% 



Table 1: Estimated contributions to the tagged and untagged samples for the three different 
charmonium decays and the non-resonant backgrounds (see text). 

mass window of ±150 MeV around the nominal charmonium masses is determined from the 
side-band events. It is found to be negligible for the ^^(25*) cascade decay samples. The elastic 
and proton dissociative contributions in Table [l| are determined using a Monte Carlo simu- 
lation of the forward region of the detector and the beamline. In the simulation the relative 
contributions from the proton dissociative and elastic channels are adjusted to reproduce the 
distributions measured in the forward detectors. The uncertainty on the determination of the 
admixtures is evaluated by variation of the simulated forward detector efficiencies and of the 
simulated cross section dependence on the mass My of the dissociated proton system p5|]. The 
resulting uncertainty on the final cross sections is 6%. Within experimental errors, equal cross 
sections for elastic and proton dissociative photoproduction of charmonium are found for the 
measured My and t regions, which is in agreement with previous results [[I|,pl,p^]. 



The fraction of erroneously tagged events in the elastic sample increases with \t\ and, with 
the present analysis method, no distinction between elastic and proton dissociative events can 
be made for |t| > 1 GeV^. However, since the t-dependence of the elastic channel is much 
steeper than that for proton dissociation, the elastic contribution at > 1 GeV^ can safely be 
neglected. 



2.5 i/^{2S) and J/i/j Signal Determination 

For the measurement of the W^p dependence the numbers of signal events are determined for the 
diffractive samples in four bins of Wyp. For events with direct decays into muons a simultaneous 
fit is used in which the ^(25) or J/tp signals are parametrised by Gaussian distributions and 
the non-resonant background follows a linear dependence. In an alternative method the non- 



resonant background is estimated using the generator LPAIR [ |20| ] and is subtracted from the 
number of signal events, counted in mass windows of ±150 MeV width around the nominal 
J/-ip and i'{2S) masses. For the number of J/-ip signal events agreement within 2% between 
the two methods is found. For the direct ipi^S) decays into muons the number of signal events 
obtained from the fits differs by up to 10% from that obtained when using LPAIR to describe 
the background. These differences are taken as systematic errors. 

For the direct decays of '0(25') and J/ip into electrons the number of signal events is determined 
using a fit to the mass spectrum with Gaussians for the J/ip and ip{2S) signals, an exponential 
distribution for the radiative tails and a linear dependence for the non-resonant background. 
The systematic error in this method is estimated by comparing the result for the non-resonant 



background with that of the Monte Carlo generator GRAPE \^A\\. For J/^ decays into electrons 
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the resulting systematic error is 5%. For the direct 4'{2S) decays into electrons the uncertainty 
on the signal determination is estimated to be 15%, which is larger than that for J/ip due to the 
larger background fraction. 

For the cascade decays the non-resonant background is very small and all events found in the 
range |M«^7r — Mu — 0.59 GeV| < 0.06 GeV are assigned to the signal (see Fig. ^jc). The 
number of J/ 4' events is corrected for the fraction of ^{2S) events with decays into a J/'ip and 
neutral particles. This correction is estimated to be (3.4 ± 0.9)% based on the results of this 
analysis (see below). The contribution to the ^'(25') sample from decays of ^{•iS) and higher 
excited states is expected to be small and is neglected. 



3 Results 

The results are given for an average (Q^) ~ 0.055 GeV^ and cover an energy range 40 < W^p < 
150 GeV. For the proton dissociative channel the kinematic region is restricted to \t\ < 5 GeV^ 
and (My/W^p)^ < 0.05. 



3.1 Energy Dependence of cr{il;{2S))/cr{J/'il;) 

The ratio of ^(25*) to J/4' cross sections R = a{tl){2S)) / o{J / ij)) is measured as a function of 
W^p. The total diffractive ^'(25') and J /ip samples, i.e. the sum of the tagged and the untagged 
samples as defined in section |2|^, are used in order to minimise the systematic and statistical 
errors. The experimental signatures of the leptonic decay channels of ^(25*) and J /i}) mesons 
are very similar and the large sample of J /i}) mesons has been used to support the study of the 
experimental systematics for the ip{2S) measurement. In the context of this analysis a complete 
measurement of the J /i}) cross section was performed [pS]], yielding a W^p dependence which 
is consistent with the results of our earlier publication [^. 

The cross section ratio is obtained by taking the ratio of the corrected numbers of ^(25") events 
to the corrected number of J /ip events in each ^(25) decay channel separately. Each of the 
samples is corrected for its specific trigger efficiency, event selection efficiency and the ^^(25*) 
and J branching ratios[|. In the evaluation of the systematic error on the cross section ratio the 
errors on the integrated luminosity, the trigger efficiency, the lepton identification and track re- 
construction efficiencies and the detector acceptance are found to be highly correlated between 
the J/ip and the il){2S) samples and therefore largely cancel. The remaining errors from these 
sources amount to 5% in total. For the direct decays the uncertainty in the determination of the 
number of V^(2S') signal events is the dominant systematic error contributing 10% for decays to 
/i+yu" and 15% for decays to e^e" (see section ^31). For the cascade decays the uncertainty in 



the pion track reconstruction efficiency of 12% dominates the systematic error. 

Figure |^ and Table || show the measured cross section ratio as a function of W^p. In the figure 
the inner error bars reflect the statistical errors, the outer error bars show the statistical and the 



^The branching ratios for the V'(25) are (0.79 ± 0.05)%, (0.77 ± 0.17)% and (34.8 ± 2.8)% for the decays to 
e+e~, /i^'^/i^ and J/i/'tt+tt^, respectively. For the J/?/' decays into muons and electrons the branching ratios are 
(5.88 ± 0.10)% and (5.93 ± 0.10)%, respectively M. 
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Figure 5: The ratio R{W^p) = (y^(2S)/(^j/ti} for events with z > 0.95. The inner error bars show 
the statistical error. The outer error bars show the statistical and systematic errors added in 
quadrature. An additional normalisation uncertainty of 0.007 due to the errors on the branching 
ratios is not shown. A fit i? oc {Wyp/90 GeY)'^^ with a value for A6 of 0.24 (solid line) and 
predictions from (dashed-dotted line) and [ |T0| ] (dashed line) are also shown. 



systematic errors added in quadrature, excluding the uncertainty from the branching ratios. The 
overall ratio, shown in Fig. ^ as a dotted line, yields the value 

R = 0.166 ± 0.007(stat.) ± 0.008(s?/s.) ± 0.007(5i?), 



in good agreement with an earlier measurement [|T2|]. This result is obtained by averaging over 
the ratios for the different decay channels and integrating over the four bins in W^p, taking 
proper account of the correlations in the branching ratio errors. The solid line shows a fit of the 
form R{W^p) oc W^p^ to the data points yielding a value of A5 = 0.24 ± 0.17 with = 3.1 
for 2 degrees of freedom, where the error includes systematic and statistical uncertainties added 
in quadrature. The energy dependence of the diffractive ^(25*) photoproduction cross section 
is thus similar or possibly slightly steeper than that for J/-ip mesons. The calculations of ^ 
and of [[l^]Q which are valid for elastic photoproduction are also shown. Both groups predict a 
somewhat steeper energy dependence for the ^(25) than for the J/^p meson, compatible with 
the data in both slope and normalisation. 



3.2 t-Dependences of the Elastic and Proton Dissociative Channels 

In order to study the t-dependences, the data are divided into the tagged and untagged samples 
as described in section Figure ^ shows the normalised differential cross sections 1/a-da/dt 
for events in mass windows of ± 150 MeV around the nominal J/ip and i'{2S) masses for each 

"^The different parametrisations given in [ [l0| ] vary considerably in normalisation. In Fig. ||the parametrisation 
referred to in [0 as 'GBW(BT)' is shown. 
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Figure 6: Normalised differential cross sections 1/a ■ da/dt as a function of |t| for events a) 
in the J/ 4' mass window without a signal in the forward detectors (untagged), b) in the J/ib 
mass window with a signal in the forward detectors (tagged), c) untagged and d) tagged events 
in the 11^(23) mass window for direct decays to lepton pairs, e) untagged and f) tagged ip{2S) 
events with cascade decays. The solid lines show the results of the fits described in the text. The 
dashed (dotted) curves show the contributions from the elastic (proton dissociative) processes, 
respectively. For the direct decays into leptons (a-d) the contributions from the non-resonant 
background (dashed-dotted curves) are also shown. The shaded bands indicate the fit regions. 
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of the three samples J /ip, ip{2S) with direct and 4'{2S) with cascade decays. The data samples 
consist of events with decays into both muons and electrons and are corrected for efficiencies 
but not for remaining backgrounds. The untagged (tagged) cross sections are displayed in the 
left (right) column of the figure. The main contribution to the untagged cross section is from the 
elastic diffractive process with contaminations from the untagged proton dissociative process 
and from non-resonant background. The main contribution to the tagged cross section is from 
the proton dissociative process, with contaminations from the non-resonant background and 
from the elastic process. 

The method used to analyse the t-dependence is optimised for samples with large non-resonant 
background, i.e. the direct 'ip{2S) decays. The same method is used for the analysis of the 
t-dependences of the J/tp samples and the ^(25) samples with cascade decays, for consis- 
tency. The t-dependences of the elastic and the proton dissociative contributions are described 
by single exponential distributions e''*. The slope parameters b are extracted using an iterative 
procedure of combined fits to the tagged and untagged samples. The non-resonant backgrounds 
are described by sums of two exponential functions. The t-dependence and relative normalisa- 
tion of the non-resonant background are determined from an analysis of the side-band events 
in the range A < Mu < 8 GeV, and are fixed in the fit procedure. The observed behaviour is 
reproduced by the LPAIR and GRAPE generators. 

In the iteration procedure the relative fractions of the background contributions are fixed to the 
numbers given in section ^!4| . In the first step a fit is performed for the proton dissociative 
slope parameter in the tagged sample, neglecting the elastic contamination. This estimate is 
then used to compute the proton dissociative contamination in the untagged sample. The Monte 
Carlo simulation is used in this calculation to account for the t-dependence of the probability 
of failing to tag proton dissociation events. Secondly, a fit to the untagged sample is performed 
to extract the slope parameter of the elastic channel, correcting for the untagged proton disso- 
ciative contribution. This measurement of the elastic channel is then used to estimate the elastic 
contribution in the tagged sample and to improve the measurement of the proton dissociative 
process obtained in the first step. The fit procedure converges after one iteration since the elas- 
tic contamination in the tagged sample is small (see Fig. The elastic slope parameters bei 
are extracted for each of the three untagged samples in the region of 0.07 < \t\ < 0.9 GeV^ 
excluding the regions with large non-resonant background (lowest \t\ bin) or large contribution 
from proton dissociation (highest \t\ bin). Similarly, the proton dissociative slope parameters 
bpd are extracted from the tagged samples in the region of 0.15 < \t\ < 3 GeV^. The lowest 
\t\ bin is excluded since for proton dissociation, the minimum kinematically allowed value of 
\t\ can be large for large masses My. The highest \t\ bin is excluded because at large values of 
|t| > 3 GeV^, the data deviate from an exponential behaviour. In Fig. ^the solid lines indicate 
the parametrisations of the data as obtained from the fits, including the contributions from the 
signals and the backgrounds. 

The systematic uncertainty on the b values is estimated by varying the selection cuts and trigger 
conditions, the |t| -dependences of the efficiencies and acceptances and the admixtures of the 
background contributions (both relative normalisations and shapes) resulting in the uncertain- 
ties given in Table |3|. The dominant error source for the elastic slope parameter is the proton 
dissociative background. For the proton dissociative slope parameter the t-dependences of the 
trigger and event selection efficiencies are the dominant error sources. 
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From the untagged samples the fits yield the elastic slope parameters b^i"^^^ = (4.31 ± 0.57 ± 
0.46) GeV-2 for the ij{2S) and b^J^ = (4.99 ± 0.13 ± 0.39) GeV^^ for the J/^, where the 
first error is statistical and the second error describes the systematic uncertainties. The result for 
b"^!^ confirms the previous measurement [|^]. From the tagged samples the proton dissociative 
parameters fej'f^^ = (0.59±0.13±0.12) GeV-^ for the ^(25) andb-^J/ = (1.07±0.03±0.11) 
GeV^^ for the J/ip are obtained. These results are listed in Table ^ together with the results for 
separate fits to the samples with muons or with electrons alone and to the samples with direct 
and cascade ^(25*) decays. Good agreement is observed between the results from the different 
samples. 

The measured slope parameters bpd with proton dissociation are considerably smaller than in 



the elastic case. This is in agreement with expectations [|15|]. The t-dependence of elastic 
ip{2S) photoproduction is similar to that of J/^j- This is in agreement with expectations in the 
additive quark model [ p31 ] and also within the colour dipole model [Q. For proton dissociation a 
somewhat shallower t-dependence is measured for the ip{'2,S) than for the J /ip. The difference 
between the proton dissociative slope parameters of J /if) and ^'(25') mesons amounts to 2.3 
standard deviations and can be interpreted as an effect of the ijj{2S) wavefunction node on the 
t-dependence. In proton dissociation this is more visible than in the elastic process due to the 
lower overall values of the slope parameter bpd- 



4 Summary 

The W-yp and t-dependences of diffractive i'{2S) photoproduction have been measured for the 
first time. For the total diffractive cross section ratio R = cr^{2S)/<^j/ip ^ value of 

R = 0.166 ± 0.007(stat.) ± 0.008(s?/s.) ± 0.007(5i?) 

is obtained in the region 40 < VF^p < 150 GeV and < 1 GeV^ consistent with and improving 
the errors of our previous measurement [|T^. The cross section ratio has been measured in four 
bins of Wyp and the result indicates that the energy dependence of tp{2S) photoproduction is 
similar or possibly somewhat steeper than that of J/ip production. The data are well described 
by pQCD calculations 



The t-dependences of elastic and proton dissociative charmonium production have been mea- 
sured. The t-dependence of elastic ^(25*) meson production is compatible with an exponential 
behaviour oc e''* in the measured range of |t| < 0.9 GeV^, with fej/^^^ = (4.31 ± 0.57 ± 0.46) 
GeV-2. This value is similar to the resuk for J/ip mesons, 6^/*^ = (4.99 ± 0.13 ± 0.39) GeV'^ 
and thus confirms the predictions of [@,|l3]- The proton dissociative slope parameter has been 
determined to be b^^^^^ = (0.59 ± 0.13 ± 0.12) GeV^^ for the 4j{2S) meson, somewhat smaller 
than that for the J/ip meson, which is measured to be fo^^^ = (1.07 ± 0.03 ± 0.11) GeV"^ 
These results are well described by pQCD calculations which take into account the differences 
in the wavefunctions of the J/ ijj and its radial excitation. 
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W^p Interval [GeV] 


(Wyp) [GeV] 




40 - 70 


53.2 


0.150 ±0.013 ±0.011 


70-90 


79.4 


0.181 ±0.015 ±0.013 


90 - 110 


99.5 


0.149 ±0.015 ±0.012 


110-150 


128.3 


0.203 ±0.018 ±0.014 



Table 2: Results for the ip{2S) to J/ip photoproduction cross section ratio R — cr^{2S)/<^j/ip in 
four bins of Wjp together with the statistical and systematic errors. An additional normalisation 
uncertainty of 0.007 due to the branching ratios is not included in the errors. 









i^{2S) - 


-^££7171 


■J/i'- 




Systematic errors [GcV^^] 


Abel 


Abpd 


Abel 


Abpd 


Abel 


Abpd 


Trigger and selection efficiencies 


0.2 


0.1 


0.2 


0.1 


0.2 


0.1 


^-dependence of forward tagging 


0.2 


0.03 


0.2 


0.03 


0.2 


0.03 


Normalisation of non-res. background 


0.11 


0.01 










t-dependence of non-res. background 


0.05 


0.04 










Normalisation of p-diss.(elas.) background 


0.35 


0.01 


0.35 


0.01 


0.27 


0.01 


t-dependence of p-diss.(elas.) background 


0.06 


0.05 


0.06 


0.05 


0.03 


0.03 


Total error 


0.47 


0.12 


0.45 


0.12 


0.39 


0.11 



Table 3: Table of systematic errors (in GeV~^) in the determination of the elastic and proton 
dissociative slope parameters h of '^{'2,3) production with direct and cascade decays and J/V' 
production. 



ip{2S) — > ii/i 
^{23) ee 
jp{23) fxfxn+n- 
ip{23) — > ee7r'*"7r~ 
J/ip ^ fifi 
J/i{j ^ ee 



bei[GeV- 
4.76 ± 0.78 
3.51 ± 2.44 
3.19 ±0.96 

5.91 ±2.74 
4.98 ±0.15 
5.05 ±0.25 



0.69 ±0.22 
0.42 ± 0.35 
0.53 ±0.28 
0.57 ±0.26 
1.10 ±0.04 
1.01 ±0.04 



ij{23) 
i'{23) 



TT 



4.69 ± 0.73 ± 0.45 
3.88 ± 0.92 ± 0.47 



0.62 ±0.18 ±0.12 
0.56 ±0.19 ±0.12 



4.31 ±0.57 ±0.46 
4.99 ±0.13 ±0.39 



0.59 ±0.13 ±0.12 
1.07 ±0.03 ±0.11 



Table 4: Slope parameters for charmonium production for the elastic (left column) and pro- 
ton dissociative (right column) component. For the separate electron and muon channels only 
the statistical errors are given. For the results combining both channels, both statistical and 
systematic errors are quoted. 
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